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Abstract 

Suppose that m drivers each choose a preferred parking space in a linear car 
park with n spaces. Each driver goes to the chosen space and parks there if it is 
free, and otherwise takes the first available space with a larger number (if any). If 
all drivers park successfully, the sequence of choices is called a parking function. In 
general, if k drivers fail to park, we have a defective parking function of defect k. 
Let cp(n, m, k) be the number of such functions. 

In this paper, we establish a recurrence relation for the numbers cp(n, m, k), and 
express this as an equation for a three-variable generating function. We solve this 
equation using the kernel method, and extract the coefficients explicitly: it turns 
out that the cumulative totals are partial sums in Abel's binomial identity. Finally, 
we compute the asymptotics of cp(n, m, k). In particular, for the case m = n, if 
choices are made independently at random, the limiting distribution of the defect 
(the number of drivers who fail to park), scaled by the square root of n, is the 
Rayleigh distribution. On the other hand, in the case m = u>(n), the probability 
that all spaces are occupied tends asymptotically to one. 



1 Introduction 

A car park consists of n numbered spaces in a line. The drivers of m cars have inde- 
pendently chosen their favorite parking spaces. Each driver arrives at the car park and 
proceeds to his chosen space, parking there if it is free. If the chosen space is occupied, 
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the driver continues on towards the larger-numbered spaces and takes the first available 
space if any; if no such space is available, the driver leaves the car park and goes home. 
What is the probability that everybody parks successfully? Equivalently, how many of 
the n m sequences of choices by the drivers lead to everyone parking? (Such a sequence is 
called a parking function.) 

This problem was first raised in the 1960s in connection with hashing [11]. In the case 
of m = n, a short and elegant proof of the formula (n + l)™" 1 was given by Pollak [7]. 
From these beginnings, a substantial theory of parking functions has been developed, 
with links to trees (as one would expect from the formula above) and priority queues [9], 
partitions [20], and representation theory [10]. More recently, generalizations of parking 
functions have found application in areas like the modelling of percolation [14] , the Abelian 
sandpile model [17] and branching processes [5]. 

In this paper, we are concerned with the probability that k drivers fail to park suc- 
cessfully. We call the corresponding assignments a defective parking function of defect k. 
Suppose that m cars attempt to park in a linear car park with n spaces according to the 
above rules; let cp(n, m, k) be the number of choices which result in exactly k drivers 
failing to park. 

The concept is related to that of x-parking functions introduced by Stanley and Pit- 
man: For a tuple of integers x = (x±, . . . ,x n ) with n G N, an x-parking function is a 
sequence (ax, ... , a n ) whose ordered permutation (pi, ... , b n ) satisfies hi < x\ + . . . + x.- t 
for all 1 < i < n. 

These generalized parking functions have been extensively studied. Stanley and Pit- 
man related their number to the volume polynomials of certain types of polytopes and to 
that of plane partitions [16]. Later, Yan and Kung investigated moments of sums of their 
numbers by using Goncarov polynomials [12, 13]. 

Classical parking functions correspond to (1,1,..., l)-parking functions while defec- 
tive parking functions of defect k correspond to (n—(m — k) + 1,1, ... ,1,0, ... , 0)-parking 
functions which in addition are not (n — (m — k), 1, . . . , 1, 0, . . . , 0)-parking functions. For 
tuples x = (a, b, . . . , b, c, 0, . . . , 0), an explicit formula was derived in [16] using the con- 
nection to polytope volumes and was reproven by Yan with combinatorial means in [21]. 

We establish a new recurrence relation for the number of defective parking functions, 
allowing us to formulate an equation defining the corresponding three- variable generating 
function. Applying the kernel method [6, 18], we solve this equation explicitly, and then 
extract the coefficients. We reobtain the fact that the cumulative totals turn out to be 
partial sums in Abel's binomial formula [1] as shown in [16] and [21] within the context of 
x-parking functions. In fact, the parking function approach may be used to prove special 
cases of this identity. 

We then investigate the asymptotical behavior of defective parking for prominent cases. 
Spencer and Yan [19] have studied asymptotics of parking functions with a defect equal 
to the difference between the number of cars and the number of spaces (the case in which 
all parking spaces end up being taken). We extend these results to the case of arbitrary 
defects. In particular, we include the case in which the number of cars is less than the 
number of spaces. 
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First, we show that, for fixed k and £, the limit of cp(n, n + £,k)/ cp(n, n + £, 0) exists, 
and compute its value. For example, the limiting value of cp(n, n, 1)/ cp(n, n, 0) is 2e — 3. 

To survey the limiting shape of the distribution, we need appropriate scaling, which 
turns out to be by the square root of n. We show that, if m = n + \_y\fn\, then the 
limiting probability of at most [^V^J drivers failing to park is 



a surprisingly simple result, given the complicated form of the exact formula. 

For y = (that is, in the case of m — n), this limiting distribution is the Rayleigh 
distribution with parameter 1/2. (This occurs as the distribution of the length of a 
random vector in the plane whose coordinates are independent normal variables with 
standard deviation 1/2. We do not know of a direct connection of this with our problem.) 

We also investigate the limiting probability that all parking spaces are occupied. Ob- 
viously, for m strictly smaller than n this probability is zero. We show that for m = [An J 
with fixed A G M + and k = m — n, 



An alternative interpretation of the above car park problem involves a variation on the 
coupon collector problem. In the original problem, there are n distinct items. If a collector 
acquires random items, she will have approximately n/e duplicates after collecting the 
first n items and will need to collect about n ■ logn items before she has a complete set. 
But suppose the items are of strictly decreasing value and she has the option of trading 
duplicate items; each item may be traded for any other one of lower value. Ideally, she 
trades duplicates against the next most valuable item she does not yet possess. Then, 
we show that she will receive only about ^/n duplicates among the first n items and a 
complete collection already with n ■ f(n) items for any function /: N — > N for which 
lim^oo f{n) = oo holds. 

We conclude this introduction with Pollak's lovely proof, adapted to the general case: 
for m cars and n spaces (with m < n), the number of ways in which every driver parks 
successfully is (n + 1 — m)(n + To see this, consider a circular car park with n+ 1 

spaces, for which the same rules apply. Now everyone will park successfully and there 
will be n + 1 — m empty spaces; such a choice will be a parking function (for the original 
problem), if and only if space number n + 1 is empty. By symmetry, this will happen in 
a fraction (n + 1 — m)/(n + 1) of the total number (n + l) m of choices. We will see later 
that our argument reproduces this result as an essential step in the working-out of the 
kernel method. 
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Table 1: The table shows the numbers cp(n, n, k) for n — 1, . . . , 10 and k — 0, ... ,9 which 
count all car parking assignments of n cars to n spaces, such that k cars are not parked. 

2 A Functional Equation 

Let cp(n, m, k) be the number of assignments of m drivers to a car park with n spaces, that 
result in exactly k drivers leaving in the end, where the parking strategy of the drivers is 
as described in the introduction above. There are n m such assignments. Then cp(n, m, k) 
can be concisely expressed as the number of functions /: {1, . . . , m} — > {1, . . . , n} for 
which for all i £ {1, . . . , n} the set f~ l ({n + l — n}) has a size of at most k + i, and 
at least one of these sets has a size of exactly k + i. 

Values of cp(n, n, k), the case in which the number of drivers and the number of spaces 
coincide, can be found in Table 1. 

We will now derive a recursion formula by transforming the parameters so that they 
are more suitable for our purpose. Let r be the number of spaces that end up unoccupied, 
and let s be the number of occupied spaces. 

Definition 1. For r,s,k £ No let a(r,s,k) denote the number of choices for which r 
spaces remain unoccupied, s spaces are occupied in the end, and k people drive home. 

This obviously means that there are n = r + s spaces in total, and that m = k + s 
drivers arrive. Observe that cp(n, m, k) = a(n — m + k,m — k, k) is the number of 
assignments for the car parking problem with n parking spaces, m visitors, and k drivers 
going home. Correspondingly, cp(n, n, k) = a(k,n — k, k). Thus, finding a solution for 
a(r, s, k) will yield a solution for the original problem. We extend this definition to all 
integers by setting a(r, s, k) =0 whenever r, s, or k is smaller than 0. For the newly 
introduced numbers, we get the following recursive formula: 

Lemma 2. For r, s, k £ No, the number of assignments of s + k drivers to r + s spaces, 
such that r spaces remain empty, s spaces are occupied, and k drivers leave, is recursively 
defined by 

{1 if r = s = k = 0, 
a{r - 1, s, 0) + J2i=i ' a ( r ' s ~ l ) k = and (r > or s > 0), 

Zi+U^-^s-hi) i fk>0. 
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Proof. If r = s = k = 0, there exists exactly one assignment. Next, let k > 0, as in the 
third case. Since k > drivers leave in the end, there are at least k + 1 drivers that arrive 
at parking space number r + s, counting the ones that actually chose it, as well as the 
ones that did not. An assignment of the s + k drivers to the r + s parking spaces satisfies 
this condition, if and only if for some i G {1, . . . , k + 1} there are k + 1 — i drivers that 
actually choose the last space and % drivers that arrive at the space, although they have 
not chosen it. For different values of i, the corresponding assignments of the drivers must 
differ. 

There are ( fc ^^ J ways to choose the k + 1 — i drivers that actually pick the last space 
and a(r, s — assignments of the remaining s — 1 + % drivers to the first r + s — 1 
spaces, such that exactly % of them will arrive at the last space. Since these assignments 
are independent of the choice of k + 1 — i drivers that choose the last space, the claimed 
recursion holds. 

Finally, let k — but r > or s > 0. If some driver arrives at the last parking 
space, whether actually choosing it or just taking it due to lack of other spaces, the same 
recursion as for k > holds. Otherwise, the last space will be empty and the number 
of assignments in which this happens is equal to the number of ways all s drivers can be 
assigned to the first r + s — 1 spaces, such that no driver has to leave. For this, there are 
exactly a(r — 1, s, 0) ways. □ 



We restate the recursion formula from the previous lemma as 

™ ( s+k \ 

a(r, s, k) = l{ r=s=k =o}{r, s, k) + l {k=0} (k) ■ a(r - 1, s, 0) + y k + 1 _ .J 



■ a(r, s — 1, i) 



where l{k=o}(k) and l{ r=s =k=o}( r , s j k) are the characteristic functions which are one, if 
k — respectively r = s = k = and zero, otherwise. Dividing both sides by ( s ~\ t k ) yields 

a(r,s,k) a(r-l,s,0) s ^ (k + 1\ a(r, s - 1, i) 
^s+kj = 1 {r=s=k=o}{r, s, k) + l {k=0} {k) \ i J pi+^j — • 

(1) 

The previous equation suggests to represent a(r, s, k) by a generating function which is 
ordinary in u and exponential in a combination of v and t. 

Lemma 3. Let A be the formal power series in the three variables u, w, and t defined by 



A(u,v,t) := a(r,s,k) 



V 3 t K 



(s+fc)! • 
r,s,fc>0 



Then A is the unique solution of 



= (2 e * - 1) • A{u, v, t) + (u - f ) • A(u, v,0) + l 



in the ring of formal power series in u, v and t. 
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Proof. Multiplying both sides of Equation (1) by u r ^y and summing both sides over 
the parameters r, s, and k, followed by the usual manipulation, such as index-shifts and 
factorizing the product of A(u, v, t) and e', shows that the definition of A and the equation 
stated in this lemma are equivalent. □ 



3 An Explicit Formula 

We will now proceed to find an explicit formula for the coefficients of the generating 
function. In general, equations like the equation in Lemma 3 cannot be solved directly, 
since both A(u, v, t) and A(u, v, 0) are unknown. Instead, we resort to using the so-called 
kernel method [6, 18]. Writing the equation in Lemma 3 as 

K(v, t) ■ A(u, v, t) = (u - f ) ■ A(u, v, 0) + 1 

with the kernel K(v, t) = 1 — |e 4 , we solve for A(u, v, 0) by setting the kernel equal to zero, 
which is here equivalent to finding a formal power series t(v) for which K(v,t(v)) = 0. 
The solution to t = ve l is the well-known tree function t = T(v) enumerating rooted trees 
on % labelled nodes, which is standardly expressed in terms of the Lambert W-function [3] 
as T(v) = — W(— v ) and has series expansion 

°° -i-i 
i=i 

Therefore 

e T(v) 

A(u,v,0) = T —^, (2) 

which can be substituted into the equation of Lemma 3 to derive an explicit expression 
for A(u, v,t). This proves the following lemma. 

Lemma 4. The generating function for the car parking problem is given by 

1 u — j e T ^ 

Applying Lagrange inversion to Equation (2), we obtain the explicit expression 

A(u, v, 0) = ( r + *) • ( r + s + • uTV i ■ ( 3 ) 

r,s>0 

The coefficients of A(u, v, 0) have already been obtained in the introduction by a direct 
combinatorial method. It becomes apparent that we can express the car parking numbers 
cp(n, m, k) = a(n — m + k,m — k, k) in terms of the following sums, which also have a 
direct combinatorial interpretation. 
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Definition 5. For n,m,k G No, let 
S(n, m, k) := 



n m if k < m — n, 

m—k 

y"] ('") • (n—m+k) ■ (n—m+k+i) 1 ' 1 ■ (m—k—i) m ~ l otherwise. 

i=0 



The car parking numbers then calculate as follows. 

Theorem 6. Let n,m,k G No- Then, the sum S(n, m, k) counts the number of car parking 
assignments of m cars on n spaces, such that at least k cars do not find a parking space, 
that is, 

m 

S(n, m, k) = cp(n, m, j) . 

j=k 

Equivalently, the car parking numbers cp(n,m,k) are given by 

cp(n, m, k) = S(n, m, k) — S(n, m, k + 1) . 

Proof. Expanding the explicit form of A(u, v, t) in Lemma 4 with the help of Equation (3) 
leads, after some lengthy manipulations, to 

s 

s 3+k v s f.k | ^ > \ \ \ fs+k\ _ „ _ i 1 . / j\s+k—i _ v s t k 



A(«,.,t)=EEw+EEEE(T*) •'■•('■+•■) 



(s+fc)! 

s>0 k>0 r>l s>0 fc>0 i=0 



r>0 s>l fc>0 i=0 



+k—i _ r v B t k 



v s t K 



= {S{r + s,s + k,k)-S(r + s,s + k,k+l)) 

r,s,k>0 

From this, we read off directly that 

a(r, s, k) = S(r + s, s + k, k) — S(r + s, s + k, k + 1) . 

The statement of the theorem follows from the relation between cp(n, m, k) and a(r, s, k). 

□ 

Note that we can rewrite S(n, m, k) for k > m — n as 



m—k / \ 

S(n,m,k) = ( • ) C P 



(n — m + k + % — 1, i, 0) • (m — k — i) r 



This observation [8] leads to a direct combinatorial proof of Theorem 6. 
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Alternative proof. If at least k cars do not find a parking space, then there are at least 
£ = n + k — m empty parking spaces. Assuming the £-th empty space occurs at position 
£ + i, there are i cars successfully parked in the £ + i — 1 spaces to the left of it, which is 
counted by cp(£ + % — 1, i, 0). Selecting these i cars out of all m cars can be done in ( m ) 
different ways. The remaining m — i cars are assigned to the n — i — i rightmost spaces 
in in — £ — z) m_ * different ways. Summing over all possible values of i leads to 

n-i , s 

S(n, m, k) = ( m ) ■ C P(^ + i-l,i,0)-(n-£- i) m_i . 

□ 

Another way to derive this theorem is to reinterpret defective parking functions in 
terms of x-parking functions as described in the introduction and then apply the results 
from [16]. 



4 Abel's Binomial Identity 

In the last section we saw that the number of assignments of m cars to n spaces, such 
that at least k drivers fail to park, is the sum S(n,m,k). Interestingly, S(n,m,k) turns 
out to be a partial Abel-type sum as it appears in Abel's Binomial identity [1]. 

Lemma 7 (Abel's Binomial Identity). 

m / \ 

( m ) ■ a ■ (« + <_1 • ( b - m_i = (a + b) m for alla,beR,men 



ci- 



In fact, our approach gives a proof of this identity for the case a, b, m e No and b — m. 
(Put k = 0, a = n — m and b = m in the defining equation of S(n, m, k) in Definition 5). 
We use this identity to find the following short expressions of S(n, m, k) and S(n, n, k). 



if k < m—n 



S(n,m,k) = { i-k fm 



( ■ ) ■ • (n-m+k) ■ (k-iy ■ (n+k-i)™- 1 -* otherwise, 

i=o 

5(n, n,k) =n n -J2( n ) ■ C" 1 )* ' fc " O^)* " (^+^-«) n " 
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We now give S(n, n, k) for values of k close to zero or n: 
S (71,71,0) = n n , 
S(n,n,l) =n n - (n + l) n ~ l , 

S(n,n,2) = n n -2- (n + 2) n ~ l + 2 ■ n ■ (n + l) n ~ 2 , 

S(n, n, 3) = n n - § • n ■ (n - 1) • (n + l) n ~ 3 + 6 • n ■ (n + 2) n ~ 2 - 3 • (n + 3)"" 1 , 



S(n, n,n-3) = 3 n + n ■ (n - 3) ■ 2 n ^ + | • n ■ (n - l) 2 • (n - 3) 
S(n,n,n-2) =2 n + n- (n-2) , 
S(n, n, n — 1) = 1 , 

1 ifn = 0, 
otherwise . 



>S(n, n, n) 



It follows that for m = n + £ with k > £ and £ £ Z, 

k-i 



lim n • (^lAl _ A = _ (fc _ £) . ^ Lll! . (fc _ i} * . e ^-^ =: k) 

n^oo \ ri / L — ' 2' 

v 7 i=0 

(This limit is trivially zero if k < £.) This implies that for £ < the limit 

.. cp(n,n + £,k) <f)(£,k)-(J)(£,k + l) 

lim — — r = lim 



rwoo C p(n, n + £, 0) n^oo 0(£,O) - 1) 

is finite (for £ > the denominator is zero). For example, we find 

cp(n, n, 1) , , cp(7T,,n, 2) 9 

lim , ; = 2e - 3 and lim ; ; = 3e 2 - 8e + 7/2 . 
n~*oo cp(n, rz, 0) «-»oo cp(n, n, 0) 

We conclude that in a random instance, we do not expect the number of drivers having 
to leave to be bounded by a constant. 

Abel's binomial identity also gives us a first bound on S(n,m,k). We obtain it by 
bounding the binomial coefficient appearing within the sum: Since for k > m — n, we 
know that 

m—k s j^. 

S(n,m,k) = S ^h(n,m,k) ■ I . ) ■ (n - m + k) ■ (n - m + k + i)' 1 ' 1 ■ (m - k - i) m ~ k ~ l 
i=o \ 1 J 

with 

(m — k — i) k ■ (m — k — i)\ ■ ml ml 



h(n, m, k) = — ■ — — < 



(m — i)l ■ [m — k)l (m — k)l 

It follows that 



™ m-k 



S(n,m,k) < -ttt ■ n 

[m — ky. 

We have seen how the formula for S(n,m,k), derived by Abel's identity, allows us 
to obtain asymptotic results for fixed values of k. In the following section, we analyze 
the case in which k grows as a function of n. Some of the upcoming results may be 
elementarily proven using the previous inequality. 
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5 Asymptotics 

For the parking problem with m cars and n spaces, the probability that at least k drivers 
cannot park their cars is S(n, m, k)/n m with S(n, m, k) as defined in Definition 5. The case 
k = corresponds to partially filled hash tables and has been analyzed at length in [2], 
with the most interesting asymptotic behavior obtained for n — m = O(^fri). Similarly, 
we find non-trivial behavior in the regime where both n — m and k are of order 0(^/n). 

Theorem 8. Let x G R + and y G KL Then, the limiting probability that in a random 
assignment of n + [z/a/^-J drivers to n spaces at least \ x\/n\ drivers fail to park is 

S(n,n+ \yy/n\,\xy/n\) \ e - 2x ( x -v) ifx>y, 
hm 



n 



n+lyvn] i otherwise. 



Proof. Let p(n, m, k) = S(n, m, k)/n m . First note that the case x < y corresponds to the 
case m > n + k, where p(n, m,k) = 1. The case x > y corresponds to the case m < n + k, 
where 



m—k 



p(n, m, k) = p(n, m, k, 



i=Q 



with p(n, m, k, i) given by 



p(n, m, k, i) = — ■ ( . ) ■ (m — k — i) m 1 ■ {n — m + k + i) 1 1 • {n — m + k) . 



A straightforward but slightly tedious calculation establishes that for a G [0, 1] 

t / | / — / — \ x-y ( {x-{l-a)y) 2 \ 

hm np(n, n + yy/n, Xy/n, an) = —===== ■ exp — r — 

n ^°° A/27ra 3 (l - a) \ 2a(l - a) J 



As we have uniform convergence to a bounded limiting function for a G [0, 1], it is 
permissible to approximate p(n,n + [y\/n\, [x\/n\) for large n by an integral as follows: 

n— \_x^/n\ 

lim p(n, n + [yVn\ , [x^fn\ ) = lim p(n, n + [yy/n\ , [^v^J , i) 



n— >oo 

i=0 



lim / n ■ p(n, n + y\/n, x\/n, an) da 

woo Jo 

x-y ( [x - (1- a)y) 2 \ 

exp — — da . 



^2^3(1 - a) V. 2a(l-a 

x+u(x—y) ' 



Under the substitution a = — t^Htt, this integral simplifies to 



1 f 00 / x(x-y) ( {l + u) 2 \ 

-= ■ / \ ■ exp — x ■ [x — y) ■ du = exp(— lx[x — y)) . 

/27T Jo V u A V 2u J 



□ 
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Figure 1: This figure shows a comparison between the probabilities p(n,m,k) = 
cp(n, m, k)/n m that m cars park randomly on n spaces, such that k drivers fail to park, 
and their asymptotic approximation for n = 100 and m = 90 (left), m = 100 (middle) 
and m = 110 (right), respectively. 



This theorem implies that for m < n + k, a good approximation is given by 

Cp(n,m,fc) *l.(2k-m + n)- e -^-m + n)/n 
n m n v ' 

Figure 1 shows a comparison between cp(n, m, k)/n m and this approximation for the three 
qualitatively different scenarios m < n, m = n and m > n. 

Theorem 8 also shows that in the special case m = n, where the number of drivers 
and the number of parking spaces coincide, a random assignment will result in ^Jn drivers 
leaving. In particular we get the following corollary: 

Corollary 9. Let k: N — > N. Then, the limiting probability that in a random assignment 
of n drivers to n spaces at least k(n) drivers fail to park is 

lim S ( n i n i k ( n )) _ 1° lim n ^ 0O A;(n)/- s /n = oo, 
»woo n n |l if lim^oo k{n)/y/n — 0. 

Another question one may wish to ask is how the number m of cars needs to scale with 
the number n of parking spaces to fill the car park with a finite limiting probability, and 
when this probability reaches one. Similarly, from the viewpoint of the coupon collector, 
it is reasonable to ask how many coupons are needed to obtain a complete set. 

Recall that the quantity S(n, m, m — n + 1) counts the number of car parking assign- 
ments of m cars on n spaces, such that at least m — n + 1 cars do not find a parking 
space, or, equivalently, such that at most n — 1 cars do find a parking space. Therefore, 
the probability that the car park is full is given by 

cp(n, m, m — n) S(n,m,m — n + 1) 
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We find non-trivial behavior when m depends linearly on n as stated in the following 
theorem, a similar version of which can be found in [19]. 

Theorem 10. Let A £ M + . Then, the limiting probability that in a random assignment 
of \_Xn\ drivers to n spaces all spaces are occupied is 

j. cp(n, [AnJ, \Xn\ — n) j if X < 1, 

™ - |j _ i . T(Ae~ A ) i/A>l. 

Proof. Let p(n, m) = S*(n, m, m — n + l)/?7, m . We have 

i=0 ^ ' 



2 



Choosing m = [Xn\ , we exchange summation with taking the limit and compute the limit 
term-wise, arriving at 



' A i 



lim p(n, [AnJ ) = V ^ • (i + l)*" 1 • e" A(1+i) . (4) 

n— >oo ^— ' ?! 

i=0 

Exchanging summation with taking the limit is justified by the Lebesgue dominated 
convergence theorem, as each term is bounded by the limiting expression, that is, 

< J_ . ( m ). U + . < n - 1 _ i)-- < . (i + . e -( 1 + i W" . 

The sum in Equation (4) converges for | Ae A | < 1 (in particular for all positive A). It 
can be expressed using the tree function T(v) with v = Xe~ x , and we find 

lim p(n, [Xn\) = — ■ T (Ae _A ) . 

Recalling that t = T{v) solves v = te~ l , we find that determining t = T (Ae~ A ) reduces 
to solving te~ l = Ae~ A . For A < 1, we find t = X, however, for A > 1 no such simplification 
is possible. □ 

Figure 2 shows a comparison between cp(n, m,m — n)/n m with m = [Xn\ and finite 
values of n and the limiting curve given by 1 — T(Ae~ A )/A. 

We conclude with a corollary of the previous theorem which emphasizes the threshold 
character of its statement. In this sense, the relation between the following corollary and 
the previous theorem reflects the relation between Corollary 9 and Theorem 8. 

Since e \-\_ x is an upper bound for T(Ae~ A )/A which converges to zero for A — > oo, 
almost all random assignments of m — m(n) drivers to n parking spaces will result in all 
spaces being occupied if m/n tends to infinity. 

Corollary 11. Letm: N — » N. Then, the limiting probability that in a random assignment 
of m(n) drivers to n spaces all spaces are occupied is 

cp (n, m(n), m(n) — n) JO if lim^oo m(n)/n < 1, 
n m ^ [1 if \\va. n ^ 00 m{n)/n = oo. 
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X 



Figure 2: This chart shows a comparison of the probabilities cp(n, m, m — n)/n m for 
n = 10, n = 20, and the limiting curve given by 1 — T(Ae~ A )/A from top to bottom, 
respectively. Here, m = [An J . 

6 Conclusion 

We have derived the generating function for the defective car parking numbers cp(n, m, k) 
and have used it to solve the problem of their exact and asymptotic enumeration. They are 
closely connected to the x-parking functions for which numerous structural interpretations 
exist, and these interpretations therefore transfer. However, we expect applications of our 
results in different areas, since parking functions, with or without defect, naturally capture 
various time-dependent models. 

For example, motivated by hashing with linear probing, a drop-push model for per- 
colation was proposed in [14]. Here, particles are dropped sequentially on a substrate, 
followed by the transport of the dropped particles via a pushing mechanism, caused by a 
local hard-core repulsion between particles on the substrate. If the transport is unidirec- 
tional, this is identical to the parking problem studied in this paper. 

Another example, the Abelian sandpile model, allows a decrease in the quantity of the 
system that corresponds to the number of cars. Such a decrease is not only possible but 
even necessary to prove certain stability properties of the system. Translated to parking 
functions, the notion of a defect naturally appears [4] . 

A further potential field of application is queueing theory. For instance, the branching 
process described in [5] can by viewed in this context. With respect to queueing, parking 
spaces are interpreted as time slots, at each of which exactly one task (represented by 
one car) can be processed. Whereas in a branching process, at least one unprocessed task 
must exist at every point in time; there is no reason to forbid idle time in a more general 
queueing process. 

Finally, we would like to point out that Alois Panholzer has independently obtained 
related results [15], notably a detailed description of limiting distributions. 
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